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ABSTRACT 


In  order  to  determine  the  Fourier  transform  of  a  quasi-per iodic  time 
series  (linear  problem) ,  or  the  power  spectrum  of  a  stationary  random 
time  series  (quadratic  problem) ,  it  is  desirable  that  data  be  recorded 
without  interruption  over  a  long  time  interval.  In  practice,  this  may 
not  be  possible.  The  effect  of  regular  interruption  such  as  the  day/ 
night  cycle  is  well  known.  We  here  Investigate  the  effect  of  irregular 
interruption  of  data  collection  (the  ‘breaking"  of  the  window  function) 
with  the  simplifying  assumption  that  there  is  a  uniform  probability  p 
that  each  interval  of  length  T  ,  of  the  total  interval  of  length  T  *  Nr, 
yields  no  data. 

For  the  linear  case  we  find  that  the  noise-to-signal  ratio  will 

-1  -2 

have  a  (one-sigma)  value  less  than  e  if  N  exceeds  p  (l-p)e 
For  the  quadratic  case,  the  same  requirement  is  met  by  the  less  restric¬ 
tive  requirement  that  N  exceed  p  ^"(l-p)e 

It  appears  that,  if  four  observatories  spaced  around  the  earth  were 
to  operate  for  25  days,  each  for  six  hours  a  day  (N  =  100),  and  if  the 
probability  of  cloud  cover  at  any  site  on  any  day  is  20%  (p  *  0.8), 
the  r.m.s.  noise-to-signal  ratio  is  0.25%  for  frequencies  displaced 
from  a  sharp  strong  signal  by  15  yHz.  The  noise-to-signal  ratio  drops 
off  rapidly  if  the  frequency  offset  exceeds  15  pHz. 


EXAMINATION  OF  TIME  SERIES  THROUGH  RANDOMLY  BROKEN  WINDOWS 


I.  INTRODUCTION 

In  many  astrophysics!  problems  one  is  concerned  with  the  study  of 
time  series.  It  often  happens  that  the  property  of  particular  interest 
is  the  spectrum  of  the  time  series.  In  principle,  one  may  determine  a 
time  series  to  a  prescribed  accuracy  by  making  measurements,  without 
interruption,  over  a  sufficiently  long  time  interval.  In  practice,  the 
length  of  time  over  which  the  variables  may  be  measured  will  be  limited. 
Moreover,  measurements  may  necessarily  be  interrupted  (or  otherwise  Impaired) 
for  one  reason  or  another.  The  relationship  of  the  spectrum  determined 
by  limited,  interrupted  measurements  to  the  intrinsic  spectrum  has  been 
the  subject  of  many  investigations,  as  recently  reviewed  by  Deeming  (1975). 

If  the  original  time  series  is  denoted  by  x(t),  one  may  regard  the 
measurements  y(t)  as  being  determined  by 

y(t)  »  f  (t)  x(t) ,  (1.1) 


where  f(t)  is  the  "window  function."  We  regard  x,  y  and  f  as  being  simple 
scalar  functions  but  the  procedure  may  be  generalized  to  replace  x,  y  by 
vectors  and  f  by  a  tensor. 

We  use  the  Fourier  transform  notation 


(1.2) 

(1.3) 


where  the  limits  of  integration  are  to  be  taken  to  be  - 00  to  +  00  if  other 
limits  are  not  explicitly  specified. 


If  we  are  interested  in  determining  x(&>),  the  Fourier  transform  of 


the  time  series  x(t),  then  we  may  use  the  relation 

y  (oa)  ■  Jda)’  f  )  x(<is  -  oi')  (1.4) 

to  relate  the  Fourier  transform  of  the  measured  time  series  y(t)  to  that 
of  the  original  time  series  x(t).  It  follows  from  this  experience  that  if 
f  (oj’)  has  appreciable  amplitude  at  frequencies  far  removed  from  ai’=o, 
for  example  if  it  has  a  subsidiary  peak  at  ai'-w0,  the  measured  time 
series  y (oi)  will  then  reflect  the  behavior  of  the  actual  series  not 
only  at  frequency  co  but  also  at  o>+ai0  -  The  frequencies  cu  and  thus  act 
as  aliases  of  one  another,  and  the  associated  distortion  of  the 
measured  spectrum  is  referred  to  as  aliasing. 

We  are  interested  in  the  possibility  that  f(t)  may  be  regarded  as  a 
random  variable,  expressible  as 


f(t) 


F(t;a1,o2, 


N 

,CV  =niiFn  an^  * 


(1.5) 


where  . . .a^  are  independent  random  variables  with  a  common  expec¬ 

tation  distribution.  By  the  central  limit  theorem  (Papoulis,  1965),  we 
expect  that  the  random  variable  f  (or  its  Fourier  transform  f)  will  have 
a  distribution  approximately  Gaussian  in  form  if  N  is  not  a  small  number, 
so  that  an  adequate  representation  of  f  would  be  given  by  its  mean  <f> 
and  its  standard  deviation  a(f). 
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w(oj)  »  f(o)f  (-0))  .  (1.11) 

Clearly  the  function  w(oj)  represents  the  capability  of  the  measurement 
process,  described  by  the  "window  function"  f(t),  to  determine  the  power 
spectrum  Sx(oi) .  The  function  w(<o)  may  be  expressed  in  terms  of  the  inde¬ 
pendent  random  variables 

w(co)  -  W(o>;a  ,a2,...,aN).  (1.12) 

Once  again,  unless  N  is  a  small  number,  we  expect  that  the  distribution 
of  w  will  be  approximately  Gaussian  so  that  it  may  be  characterized  by 
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its  mean  value  <w>  and  standard  deviation  cr(w). 

This  article  was  prompted  by  a  problem  related  to  the  determination 
of  normal  modes  of  oscillation  of  the  sun,  as  determined  by  measurement  of 
the  photospheric  velocity  field.  Measurements  have  been  presented  by 
Deubner  (1975)  and  by  Rhodes  et  al.  (1977),  and  their  theoretical  inter¬ 
pretation  discussed  by  Ulrich  and  Rhodes  (1977)  and  by  Ulrich  et  al.  (1978) 
For  optimum  determination  of  the  power  spectrum  of  the  velocity  field 
(expressed  as  a  function  of  wave  number),  it  is  clearly  desirable  to  make 
observations  without  interpretation  over  as  long  an  interval  as  possible. 
Away  from  polar  regions,  observations  from  a  single  station  are  interrupted 
by  the  day-night  cycle  which  leads  to  unacceptable  aliasing  of  the  data. 
Observations  made  from  a  spacecraft  in  polar  orbit  would  obviously  yield 
un-aliased  data  of  higher  quality  and  higher  frequency  resolution.  Obser¬ 
vations  made  from  the  south  pole  during  austral  midsummer  can  lead  to 
several  days  of  uninterrupted  observation  and  to  still  longer  intervals 
with  occasional,  irregular  interruption.  It  is  also  possible  to  select 
three  or  four  stations  around  the  earth  which,  in  the  absence  of  any  cloud 
cover,  could  give  continual  coverage  of  the  sun  for  many  weeks.  However, 
one  must  anticipate  that  some  of  the  data  would  be  lost  by  cloud  cover. 

It  is  clearly  desirable  that  one  should  be  able  to  make  some  esti¬ 
mate  of  the  accuracy  with  which  oscillation  modes  may  be  determined  when 
it  appears  possible  to  observe  the  sun  over  a  long  interval  of  time  losing 
some  blocks  of  time  because  of  cloud  cover.  The  purpose  of  this  article 
is  to  develop  a  model  which  enables  us  to  address  problems  of  this  type. 


After  presenting  a  few  general  formulas,  we  shall  simplify  the  problem 
considerably  by  supposing  that  observations  are  made  over  a  large  number 
N  of  equal  time  intervals,  each  of  length  x  ,  so  that  the  total  time 
interval  T  is  given  by 

T  -  NT  .  (1.13) 

With  certain  additional  simplifying  assumptions,  we  shall  consider  the 
statistical  properties  of  the  functions  f(u))  and  w(oi)  which  are  repre¬ 
sentative  of  "randomly  broken"  window  functions. 

II.  MATHEMATICAL  MODEL 

In  the  case  that  the  window  function  f(t)  is  expressible  in  the  form 
(1.5),  in  terms  of  a  number  of  random  variables,  we  wish  to  study  the  distri¬ 
bution  of  the  functions  f(cu),  w(oj),  entering  equations  (1.4)  and  (1.10). 

We  suppose  that  the  distribution  of  the  variables  a^  to  a ^  is  given  by 
the  probability  function  P(a^,  ...,  a^)  such  that  P(a^,  ...,  a^)da^  ....  daN 
is  the  probability  of  finding  in  the  range  a ^  to  a^+da^,  etc.  Then  the 
expectation  value  of  the  quantity  f(oj)  is  given  by 

<f(u»  -  PN(ct)  F(<s;a1, . . .  ,a^) ,  (2.1) 

N 

where  d  a  denotes  da.  ...,  da„,  and  P„(a)  denotes  P(a-,,  ....  avt) .  If  we 

JL  N  N  J-  N 

use  the  following  notation  for  the  variance  of  a  complex  variable  of  a 
complex  variable  z, 

a2(z)  -  o2(zr)  +  a2(Zi),  (2.2) 

where  z^  and  z^  are  the  real  and  imaginary  parts  of  z,  then  noting  that 
f(-oj)  is  the  complex  conjugate  of  f  (oj)  ,  we  see  that 
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a2  (f  (o>))  -  <f(co)  f  (-oi)>  -<f(o>»  <  f  (-  co)>  .  (2.3) 

The  first  term  on  the  right-hand  side  may  be  evaluated  from 

<w(o,)>=<f(w)  f  (-o))>»^Na  PN(a)F(u>;ct1,...>aN)  F  (- co;  (2.4) 

We  see  that  equation  (2.4)  also  gives  the  expectation  value  of  the 
"window  spectrum"  w(o>)  which  appears  in  equation  (1.10)  and  is  appropriate 
for  the  discussion  of  stationary  random  time  series.  The  variance  of  this 
function  is  given  by 


*(w)  *  <w(&)>  W  (—  6))>  -  <w(C0)>  <w(-a>)> 


(2.5) 


where 


<w(0d)  w(-  0))>  = 


/Ap„ 


(a)  ^(ojjo^,  . . 


I2  b 


*  ,CtN  |  |F(-  (i>,  , 


(2.6) 


As  indicated  in  the  introduction,  we  intend  to  consider  the  case  that 
the  observing  time  tQ  to  tN,  of  length  T,  is  divided  into  N  equal  intervals 
bounded  by  times  t^,  ...  where 


t  =  t.  +  nx 
n  0 


(2.7) 


so  that  we  may  adopt  the  form 

N 

F(t;a1,  ...,aN)  =  2  °n  h(t  ~  ‘"n  -  1}  *  h(t  ‘  tn)j  (2'8) 

n  =  1 

where  h(t)  is  the  Heavyside  function: 


h(t) 


0  if  t<0,  ) 
1  if  t  >0.  ) 


(2.9) 


We  also  assume  that  the  intervals  are  statistically  independent,  so  that 


6 


we  may  write 

PN(a)dNa  =  PCc^ . . dctN  "  jw  dal  j *  *  *  daN }  *  (2*10) 

If  we  assume  that,  for  each  Interval,  there  is  a  (uniform)  probability  p 
that  the  window  is  open  and  probability  1  -  p  that  it  is  closed,  then 

P(an)  -  p  6(an-l)  +  (l-p)S  (an)  .  (2.11) 

In  evaluating  (f(«))»  given  by  (2.1),  we  will  use 

<cxn>  -  p  .  (2.12) 

In  evaluating  the  quantity  given  by  equation  (2. A),  we  will  need  to 

a  a'  which  Is  clearly  given  by  p  if  m^  n  but  by  p  if  m  *  n. 
in  tv 

Hence 

<aa>  =  p2  +  (p  -  p2)  5  ,  (2.13) 

m  n  r  r  r  mn 

where  6  is  the  Kronecker  function.  In  evaluating  the  quantity  given 
ton 

by  (2.6),  we  need  to  evaluate  the  expectation  value  of  amanapaq*  By  con¬ 
sidering  the  various  possibilities  (m,n,p,q  all  different:  two  of  them 
the  same,  etc.)  we  find  that 

<a  a  a  a.^Wp4  +  p3-p4)  ($  +  <$  +  <5  +  &  +  6  ) 

NVn  pV  y  m  mP  m<l  n<5  pq 

+  (p2— 2p3+p4)  (6  6  +6  6  +6  5  . 

'■t'  v  '  '  mn  pq  mn  nq  mq  np) 

+  (p2-3p3+2p4)  (6npq  +  Smpq  +  Smnq  + 

+  (p-7p2+12p3-6p4)  Smpq  (2.14) 

where  6  «lifm-n»p  otherwise  0,  and  5  is  defined  similarly, 

mnp  uuipq 
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III.  EVALUATION  OF  MODEL 

For  simple  (non-random)  time  series,  equation  (1.4)  gives  the  relation¬ 
ship  between  the  Fourier  transforms  of  the  original  and  measured  time  series. 
In  this  context,  the  properties  of  the  random  window  function  f(t)  may  be 
characterized  by  ^f(co)^  and  a2( f )  . 

On  substituting  the  form  (2.8)  into  (1.5),  we  find  that 


. V  -15  £,  “n  *l»(f“T).HCo+<"'2)Tl  • 

n »  1  '  ' 

On  using  (2.1)  and  (2.12),  we  obtain 

i  1  i<o  (t  +  t  ) 

<f(<u)>  “2^  T  P  sine  (^  a>  T)  e  ^  , 

where  sine  0  *  6  ^  sin  9 . 

On  using  (2.4),  we  find  that 

<w(0))>-  (^)2  T2  p2sinc2  (|coT)  +  N_1p(l-  p) sine2 (| cot )j 
Hence,  using  (2.3),  we  obtain 


a2(f)  »  2  T2  if1  p(l  -  p)  sine2  ^  m )  . 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


For  evaluating  the  effects  of  "breaking”  of  the  window  function,  it 
is  convenient  to  normalize  the  standard  deviation  with  respect  to  the  maxi¬ 
mum  value  of  f(o>),  which  is  the  value  at  o>  *  0.  Accordingly,  we  introduce 
the  definition 


y  (a))  , 

For  the  case  under  consideration,  this  has  the  form 

S1(co)  -  N"1/2p"1/2  (l-p)1/2  (sine  |u)T|  . 


(3.5) 


(3.6) 
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For  discussion  of  the  properties  of  randomly  broken  windows  in  the 
study  of  stationary  random  time  series,  it  is  necessary  to  evaluate  the 
mean  value  and  standard  deviation  of  w(oj)  .  The  former  is  given  by  equa¬ 
tion  (3.3).  The  first  term  inside  the  brackets  has  the  same  form  as  arises 
in  the  non-random  case  (p  =  1).  The  second  term  represents  a  change  in  the 
mean  spectrum,  so  it  is  convenient  to  introduce  the  symbol  for  the 

ratio  of  the  additional  term  to  the  maximum  value  of  the  principal  term: 

A2  =  N-1  p-1  (1-p)  sine2  (-|ur)  .  (3.7) 

On  writing  equation  (2.6)  in  the  simpler  form 

<w(oj)w(- co)>  »  <F(<D;a)F(a>;a)F(-  oj;  a)F(-  <u;  a)>  (3.8) 


and  using  equations  (1.3)  and  (2.8),  we  see  that 

<w(«)w(-o»)>  =  — S  X  <Vpq>  Um*  Un*  Up  V 

(2itu)  mnpq 


(3.9) 


iojt  iojt  , 

where  U  (<o)  =  e  m  -  e  and  Um*(w)  =  U^-w)  is  the  complex  conjugate 

of  U  . 
m 

On  using  equation  (2.14),  we  see  that  this  may  be  expressed  in  the  form 

<w(w)w(-co)  >  =  — — — !p4E-  +  (p3  -  p4)  (  2^2  +  4E^) 

(  2ttco)  ^  1 


+  4  (p2  -  3p3  +  2p4) 


(3.10) 


+  (p  -  7p2  +  12p3  -  6p4)  Ec 


+  (p2  -  2  p3  +  p4)  (E6  +  2E?)j 


where 
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(That  Ej  and  E^  are  real  is  due  to  the  fact  that 
all  m.) 


2  si 


for 


Hence  equation  (3.10)  is  found  to  be  expressible  as 


<w(a>)  w  (-<o)>  -  jp4  sinc4||lj 


+  2  p^(l-p)  sine  /  wt\  sine  /  ari_ \|2  +  sine  (<uT)  +  2^1-2p)l 
N  \  2  /  \  2  sine  <oT  N  J 


♦T1  sl"'1?)[ 


p(l-p) 


/  ,  ,  sinc(ajT)  \ 
I  sine  (err)  j 


6p2  -  6p  +  1 1  | 

N2  J 


\ 

On  using  equations  (2.3),  (3.3)  and  (3.13)  and  the  definition 


(3.13) 


a  (w(co)) 


(3.14) 


we  find  that 


2  2  (to) 


2  (1-p)  sine  j  coT 
NP  '  “ 


/ ut\  sine2/  tor  \  [  1  +  sinc(a)T)  +  2(l-2p)  | 
V  2  /  \  2  /  ^  sine  (tor)  Np  J 

+(^)  -*(?)  [r 


(3.15) 
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IV.  DISCUSSION 


We  see  from  the  preceding  section  that  the  effect  of  a  random  "break¬ 
ing"  of  the  window  function  is  to  produce  an  aliasing  of  any  signal.  This 
effect  is  described  by  the  function  or  by  and  as  given  by  equa¬ 
tions  (3.6),  (3.7)  and  (3.15). 

For  the  "linear"  problem,  as  described  by  equation  (1.4)  the  mean 
Fourier  transform  of  the  window  function,  as  given  by  equation  (3.2)» 
has  the  same  form  as  it  does  in  the  non-random  case,  although  it  is 
reduced  by  a  factor  p.  The  standard  deviation  is  characterized  by  Z^, 
defined  by  equation  (3.5)  and  given  by  equation  (3.6). 

It  is  convenient  to  introduce  the  notation 

o>T  -  2T-1,  (ot=-2t"1.  (4.1) 

On  noting  that  sine  0  <  1  for  8  (  1,  and  that  |sinc  9|  <  0-1  for  0^1, 
we  see  that  Z^  <  S^  where 


The  above  equations  suggest  that  the  aliasing  is  most  severe 
within  the  range  of  a  few  times  oo  of  the  strongest  signal.  It  is 
informative  to  compare  S^(to)  with  C^(w)  defined  t>y 


<  f~(u)> 

<f~(o)> 


C1((d)  »  1,  u>  «  a) 


-  9  0)  »  GJ 

0)  "  T, 


l 

I 


(4.3) 


12 


which  represents  the  normalized  envelope  of  the  "core"  of  the  window 
function,  which  has  the  "ideal"  form  which  obtains  in  the  absence  of 
breaking.  The  variation  of  both  C^((o)  and  S^(o>)  with  03  is  shown 
schematically,  on  a  logarithmic  scale,  in  Figure  1. 

We  see  from  this  diagram  that  the  random  interruption  of  data 
collection  adds  a  "wing"  to  the  core  of  the  window  function.  This 
wing  has  a  height  ^(l-p)/Npj  over  the  range  o>c  to  a>T,  where 


Np\1/2  <*>, 


(4.4) 


provided  that  N  >  p/(l-p)  . 

We  may  readily  calculate  the  minimum  number  of  intervals 
necessary  to  ensure  that  is  below  an  assigned  level  £  for  a  given 
value  of  p.  We  see  from  (4.2)  that  S^S  £  if  N  >  where 


(4.5) 


If,  for  instance,  p  -  0.8  and  we  require  that  <  0.05,  N  must  be  at 


least  100. 


For  the  quadratic  problem  in  which  we  are  determining  the  spectrum  of 
a  stationary  random  time  series,  it  is  convenient,  as  in  our  discussion  of  the 
linear  case,  to  introduce  C2(w),  the  normalized  envelope  of  the  "core" 
of  the  window  function  for  the  case  of  uninterrupted  data  collection. 

As  we  see  from  equation  (3.3),  this  function  is  given  approximately  by 


the  formula 


C-(co)  “  1,  co<  co  , 


£)  •  - 


(4.6) 
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The  aliasing  is  now  described  by  the  functions  A2 Coj)  and  E2(<ij)  given 
by  equations  (3.7)  and  (3.15).  We  find  that  the  envelope  of  A 2 (w)  is 
bounded  by  the  function  D2(gj)  given  by 


D2(oj)  * 


AlE 

Np 


(iXdl 

T  ’ 


1=R 

Np 


(A. 7) 


The  effect  of  D2(gj)  is  to  produce  a  "wing"  to  the  window  function, 
extending  ov»t  from  oi*coc«  where  is  given  by  equation  (A«4  ) ,  as  shown 

in  Figure  3. 

Discussion  of  the  function  £2(to),  representing  the  normalized  standard 
deviation  of  the  window  function  from  the  mean  form  given  by  equation 
(3.3),  is  somewhat  more  complicated.  However,  if  Np».l,  we  find  from 
equation  (3.15)  that  E2(a>)  <  S2(<o)  where 


S2(co)  =  2(^) 

1/2 

»  0!^  0)  , 

a 

"  2i^) 

1/2 

—  .  CO<CO^U)  , 

0J  >  d 

co 

d 

.  Aze  / 

■  f 

Np  y  os 

/ 

’  T 

where 


) 

I 


(A  .9) 
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The  functions  (<w)  and  S2(co)  are  shown  plotted  in  Figure  2  for 
the  particular  case  N  *  100,  p=  0.8  and  T  *  6  hours  (a  case  discussed 
later  in  this  section).  We  see  from  this  figure  that  S 2 (w)  provides 
an  excellent  approximation  to  the  envelope  of  I2(co)  • 

The  function  S 2 (<*>)  (the  normalized  envelope  of  the  standard  deviation 
of  the  window  function)  is  also  shown  schematically  in  Figure  3.  We  see 
that  S2(o>)  ■  D2(cd)  for  co^co^.  We  also  see  that  S2(co)  <  C2(co)  for 


and  we  note  that  D„(co)  <  C  (<o)  also  in  this  range.  The  frequencies 
2  2 

0).  ,  <0  ,  ok,  are  in  the  proportion 
bed 


“b 


CO 


CO 


,-1/2 


,1/2 


(4.11) 


We  may  now  calculate  the  minimum  number  N2  of  intervals  necessary 

to  ensure  that  S2(co)  [  and  hence  E2(co)  and  A2(co)  ]  is  below  an  assigned 

level  e  for  a  given  value  of  p.  Ignoring  the  small  region  between 

co.  and  co.,  where  the  condition  is  a  little  more  restrictive,  we  see  that 
b  d 

the  required  condition  is  satisfied  if 

N  ^  .  (4.12) 

2  pe 

We  see  from  (4.5)  and  (4.12)  that  N2  is  smaller  than  by  the  factor 
G.  Hence  aliasing  is  likely  to  be  less  serious  in  the  quadratic  case  than 
it  is  in  the  linear  case. 
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JLJ-  lilLUJ, 


In  order  to  assess  the  implications  of  the  present  model  concerning 
ground-based  observations  of  solar  oscillations,  one  will  need  to  have 
detailed  estimates  of  the  expected  spectrum  (in  particular,  the  spacing 
and  relative  power  of  nearby  lines)  and  the  expected  cloud  cover  at  three 
or  four  observatories  postioned  round  the  world.  It  is  also  desirable  that 
the  present  model  should  be  extended  by  considering  separate  values  of 
p  for  each  of  the  observatories,  and  possibly  by  taking  into  account 
the  correlation  between  cloud  cover  on  consecutive  days. 

Nevertheless,  we  can  illustrate  the  results  of  this  model  by  consider¬ 
ing  a  hypothetical  situation.  Suppose  that  four  observatories  are  located 
around  the  world  in  such  a  way  as  to  give  continous  coverage  (in  the 
absence  of  cloud  cover) ,  and  that  these  observatories  are  operated  for  25 
days.  Then  N  =  100.  Suppose  that,  for  any  observatory  on  any  day,  there 
is  a  20%  probability  of  cloud  cover  so  that  p  =  0.80.  We  find  from  (4.1) 
that  (using  v  •  <i)/2tt)  ,  =  15yHz.  For  frquencies  less  than  this  value, 

(4.7)  and  (4.8)  show  that  aliasing  amounts  to  0.25%  or  less.  For 
frequencies  above  15yHz,  the  aliasing  drops  off  rapidly. 

Although  it  will  be  necessary  to  make  more  detailed  and  specific 
calculations  to  draw  definite  conclusions,  it  appears  from  the  above 
simple  example  that  it  may  be  possible  to  carry  out  high-quality  studies 
of  solar  oscillations  from  a  chain  of  ground-based  observatories. 

This  work  was  supported  in  part  by  NASA  Grant  NGL  05-020-272,  Office 
of  Naval  Research  Contract  N00014-75-0673,  and  the  Max  C.  Fleischmann 
Foundation. 
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CAPTIONS 


Figure  1.  Schematic  diagram  shoving  the  relationship  between  S^(oj)  , 

the  normalized  envelope  of  the  standard  deviation,  and  C^to)  , 
the  normalized  "core"  function,  for  the  linear  case.  The 
scales  of  both  axes  are  logarithmic. 

Figure  2.  Comparison  of  the  normalized  standard  deviation  Z2CU)  and 

the  approximate  envelope  function  S^ (00)  for  the  special  case 
discussed  in  the  text. 

Figure  3.  Schematic  diagram  showl'j,  the  relationship  between  (w)  ,  the 
normalized  envelope  of  the  standard  deviation,  02(00)  ,  the 
normalized  envelope  of  '■he  mean  displacement,  and  C^ioi)  ,  the 
normalized  core  function,  for  the  quadratic  case.  The  scales  of 
both  axes  are  logarithmic. 


Figure 


